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Abstract. Runaway electrons can be generated in tokamak plasmas if the 
accelerating force from the toroidal electric field exceeds the collisional drag force owing 
to Coulomb collisions with the background plasma. In ITER, disruptions are expected 
to generate runaway electrons mainly through knock-on collisions jT], where enough 
momentum can be transferred from existing runaways to slow electrons to transport 
the latter beyond a critical momentum, setting off an avalanche of runaway electrons. 
Since knock-on runaways are usually scattered off with a significant perpendicular 
component of the momentum with respect to the local magnetic field direction, these 
particles are highly magnetized. Consequently, the momentum dynamics require a 
full 3-D kinetic description, since these electrons are highly sensitive to the magnetic 
non-uniformity of a toroidal configuration. For this purpose, a bounce-averaged knock- 
on source term is derived. The generation of runaway electrons from the combined 
effect of Dreicer mechanism and knock-on collision process is studied with the code 
LUKE, a solver of the 3-D linearized bounce-averaged relativistic electron Fokker- 
Planck equation |2] , through the calculation of the response of the electron distribution 
function to a constant parallel electric field. The model, which has been successfully 
benchmarked against the standard Dreicer runaway theory now describes the runaway 
generation by knock-on collisions as proposed by Rosenbluth [3]. This paper shows 
that the avalanche effect can be important even in non-disruptive scenarios. Runaway 
formation through knock-on collisions is found to be strongly reduced when taking 
place off the magnetic axis, since trapped electrons can not contribute to the runaway 
electron population. Finally, the relative importance of the avalanche mechanism is 
investigated as a function of the key parameters for runaway electron formation, namely 
the plasma temperature and the electric field strength. In agreement with theoretical 
predictions, the LUKE simulations show that in low temperature and electric field the 
knock-on collisions becomes the dominant source of runaway electrons and can play a 
significant role for runaway electron generation, including in non-disruptive tokamak 
scenarios. 
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1. Introduction 

Runaway electrons have been observed in magnetic confinement fusion experiments 
during the operation of tokamaks jl]. They are also encountered in nature in solar flares 
and electric discharges associated with thunderstorms |n|. The dynamics of electrons 
in a plasma is governed by the balance between acceleration in an electric field and 
collisions with the plasma particles. Collisional friction forces acting on the electrons 
reach a global maximum at the thermal velocity {vth) and decrease for higher velocities. 
In the presence of a strong toroidal electric field {E) collisional drag may consequently 
be too weak to counteract the acceleration of electrons, which may result in continuously 
accelerated electrons, known as runaway electrons. If no other loss mechanisms than 
the collisional drag are present jS], runaway electrons may be generated if the electric 
field exceeds the critical field [7] 

HeC ^In A 
inelmoc^ ’ 

where Ue is the electron density, mo is the electron rest mass, c is the speed of light, e 
is the elementary charge, and In A is the Coulomb logarithm. The acceleration by a DC 
field of electrons that diffuse via small angle collisions beyond the critical momentum 
{pc), defined as the minimum momentum for which collisions are too weak to prevent 
acceleration of the electrons by the electric field to even higher energies, is referred to 
as the Dreicer mechanism [7]. In addition, these relativistic electrons can undergo close 
collisions with bulk electrons and transfer part of their momentum so that also the target 
electrons may get kicked into the runaway momentum region, while the momentum of 
the primary electrons remains above the critical momentum. These knock-on collisions 
can therefore lead to multiplication of the number of runaway electrons, commonly 
referred to as runaway avalanche j3]. 

Various methods to mitigate the formation of runaway electrons in tokamak plasmas 
are based on either increasing the plasma density and thereby Ec by so-called massive gas 
injection (MGI) [8], or on deconfining the runaway electrons before they can reach too 
high energy, by the means of resonant magnetic perturbations (RMP) [H]. Even though 
such mitigation methods have been demonstrated in present tokamak experiments, they 
might not provide a solution for large tokamaks like ITER uni. Therefore the formation 
of the runaway electron population is a topic in urgent need of investigation. 

Intense beams of highly energetic runaway electrons can form in tokamaks 
during plasma disruptions, fast unstable events that lead to a sudden loss of plasma 
confinement. If runaway electrons strike the first wall of the tokamak vacuum vessel the 
local energy deposition can cause significant damage [T]. Regardless of the mechanisms 
that lead to the onset of a major disruption, the post-disruption phases usually have 
similar time evolution m- They start with a fast cooling of the plasma typically 
associated with either intense radiative losses or ergodisation of the magnetic flux 
surfaces [12], referred to as the thermal quench, which occurs on a time scale on the 
order of a millisecond. Consequently the plasma resistivity p, which scales with the 
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temperature as increases rapidly. The toroidal electric field is proportional to 

the resistivity and increases dramatically in order to maintain the local current density. 
The resistive current diffusion is enhanced by the reduction of the plasma conductivity, 
such that the plasma current decays progressively. Yet, the current decay occurs over a 
much longer time scale. In this process, a fraction of the pre-disruptive plasma current 
is carried by runaway electrons. 

Disruptions are interesting but complex processes for studying the birth of runaway 
electrons, since they include magnetohydrodynamic (MHD) instabilities, anomalous 
transport and complex evolution of the magnetic field topology [13]. However, the 
generation of runaway electrons does not necessarily require the extreme conditions 
found in disruptions. In low density plasmas, the electric field can exceed the critical 
electric field also during the current flattop in a quiescent plasma, free of equilibrium 
transients, or during current ramp up or ramp down. An advantage of studying runaway 
formation in these so called non-disruptive scenarios is that the key parameters for 
the runaway electron mechanisms, mainly the electric field strength, electron density 
and temperature, can be better diagnosed than during disruptions. Runaway electrons 
have been detected in non-disruptive scenarios in several of the existing tokamaks 
na. Quiescent plasmas with nested magnetic flux surfaces are therefore more suitable 
for studying the formation of runaway electrons. In this work the formation of 
runaway electrons generated from the combined effect of Dreicer and knock-on collision 
mechanisms is studied with the code LUKE, a solver of the 3-D (one spatial and two 
momentum dimensions) linearized bounce-averaged relativistic electron Fokker-Planck 
equation |2]. The code LUKE handles arbitrary shapes of the flux surfaces, but in this 
work the magnetic flux surfaces are assumed to remain circular and concentric as in 
the Tore Supra tokamak. They are assumed to remain intact throughout the runaway 
formation process, an assumption that would be too restrictive for the thermal quench 
in disruptive scenarios. 

Modelling the evolution of the temperature and electric field in disruptions would 
require a proper description of the thermal quench including radiative or convective loss 
mechanisms and MHD instabilities. The coupling of a kinetic code capable of handling 
3-D magnetic topologies and open field lines with a fluid code such as JOREK [T^ 
would be necessary for such a purpose, but is beyond the scope of this work. The 
kinetic modelling of the formation of runaway electrons is therefore done for non- 
disruptive scenarios as found in the current flattop with constant electric field and 
plasma temperature. With the restrictions of disruption modelling in mind, the objective 
of this work is to study the formation of runaway electrons in non-disruptive scenarios 
owing to the combined effect of Dreicer and knock-on collisions with a fast solver for 
the electron distribution function, in order to make predictions for the birth of runaway 
electrons in tokamak experiments. 

The LUKE code has previously been used for current drive and Dreicer runaway 
calculations. The model uses a relativistic collision operator for small angle collisions 
and a recently added description of the large angle (knock-on) collisions leading to the 
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avalanche effect, which enables a description of the full 2-D momentum dynamics of the 
runaway population. Runaway electrons generated via knock-on collisions are typically 
scattered off with a signihcant perpendicular component of the momentum with respect 
to the local magnetic held direction. In a non-uniform magnetic held conhguration, 
highly magnetized electrons could be subject to magnetic trapping ehects resulting in 
reduced runaway electron growth rate oh the magnetic axis in comparison to estimates 
obtained for a cylindric geometry. Such toroidicity ehects are studied by implementing a 
2-D kinetic description of the knock-on momentum dynamics, including the momentum 
dynamics both perpendicular and parallel to the magnetic held lines. 

Knock-on collisions are included in the kinetic equation through a source term 
from Ref. [3], implemented along with a sink term to ensure a particle conserving 
form of the process. The bounce-averaged knock-on source term is presented in 
Sec. In Sec. the ehect of magnetic held non-uniformity is investigated. The 
role played by the magnetic mirror force on the runaway population oh the magnetic 
axis, owing to a reduction in Dreicer growth rate as well as the high magnetization 
of the knock-on electrons, is described. Finally, in Sec. the relative importance of 
the avalanche ehect compared to the Dreicer mechanism is quantihed as a function of 
plasma temperature and toroidal electric held strength. The parametric dependencies of 
the relative importance of the avalanche ehect obtained from the numerical modelling is 
related to both analytic predictions and experimental data from runaway observations in 
non-disruptive scenarios from several tokamaks. The comparison includes a low density 
hattop pulse from the Tore Supra tokamak, during which suprathermal electrons are 
observed. The analysis of this scenario supports recently published results la, showing 
that runaway electron formation requires lower density than expected from collisional 
theory, which suggests the existence of additional runaway electron loss mechanisms. 

2. Knock-on collisions model 

A knock-on collision between an existing runaway electron and a slower electron is 
considered. This paper follows the model from Ref. [3j in which the target electron is 
assumed to be at rest whereas the initial runaway travels at the speed of light in the 
direction of the magnetic held. This approximation will be justihed later in this section. 
The target electron gains a momentum p from the close collision. As both energy and 
momentum must be conserved in the collision process, the secondary electron is scattered 
with a pitch-angle with respect to the direction of the incoming electron, which cosine 
is given by the relation 



( 2 ) 


where 7 = is the relativistic factor and the momentum p is normalized to moc. 
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The relativistic electron electron differential cross section derived by M0ller |16] 
yields 


da 


1 

P7(7 - 1)^ 


5(e-r(p)), 


( 3 ) 


where r^ = e‘^/{AneomeC^) is the classical electron radius. As da/dQ decrease rapidly 
with momentum, a large fraction of secondary electrons have a moderate kinetic energy 
with 7 — 1 <C 1 and are thus scattered with a large pitch-angle characterized by 
1. Hence it is necessary to properly account for the 2-D guiding-center momentum 
dynamics in non-uniform magnetic field geometries, where the electrons are influenced 
by the magnetic trapping effect. 

The source term originally formulated in Ref. [3] is proportional to both the 
target population, i.e. the bulk electron density and the existing runaway electron 
population Ur 


S(V’,P,0 = 


Ur 


1 d 


dvrr In A p"^ dp 


1 - y/1 




( 4 ) 


where "0 is the poloidal magnetic flux surface coordinate. In the expression above, the 
collision time for relativistic electrons 
ATielmlc^ 

In A ’ 

has been introduced. 

An analytic estimate of the avalanche growth rate is obtained from integration of 
the knock-on source term in Eq. |^over the runaway region p > pc in momentum space, 
as done in the work by Rosenbluth, which yields the following expression for E > 


1 dur 1 f E \ 
Hr dt 2 rlnA / 


2.1. Implementation of knock-on collisions in the LUKE code 

The Rosenbluth model (Eq. for the runaway generation through knock-on collisions 
is implemented in the code LUKE and benchmarked against the growth rate in Eq. 
in the case of cylindrical geometry in Fig. [T| by using the same momentum thresholds 
as in Ref |3]. Electrons with a momentum larger than pre = niax[pc;p(Efc = 1 MeV)] 
are accounted for in the population n^. of primary runaways for the knock-on collision 
process. The numerical momentum grid boundary Pmax niust be chosen to be larger or 
equal to Pre and electrons leaving the domain through the boundary remain accounted 
for in Ur. 

To be valid, the Rosenbluth approximation requires that: (a) primary runaways in 
the knock-on collision process have a velocity near the speed of light, and (b) primary 
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Figure 1. The knock-on process in LUKE (crosses) benchmarked against analytic 
growth rate in Ref. [3] (dashed line), when using the same momentum space thresholds. 


electrons have a momentum much larger than target electrons. The condition (a) is 
ensured by the 1 MeV minimum condition in p^e, which corresponds to v/c > 0.94, 
whereas the condition (b) is guaranteed by restricting the model to plasmas with Te -C 1 
MeV. The Rosenbluth approximation is further justihed by the weak dependence of the 
knock-on source term upon the incident electron energy in the energy range 1 —100 MeV 

HU. 

The bulk electron density is dehned as the integral of the bulk electron distribution 
in momentum space: 



f{r,p)d^p 


ne{r). 


(7) 


The bulk and the runaway region, corresponding to p < p^e and p > Pre respectively, 
are shown in Fig. The runaway electron population is the integral over both Dreicer 
and knock-on runaway fluxes 


Urit) 



dt, 


where 7d = ff Sp(V’,P, 0 ' <78 is the integral of the particle flux through the surface 
p = Pre- In order to ensure conservation of number of particles in LUKE, a sink term is 
implemented to compensate for the knock-on source term 


S 


S+-< S+ > 


/m 

< /m >’ 


( 8 ) 


where /m is the bulk distribution, assumed to be Maxwellian and < ... >= ... d^p. 

The source and sink terms ensure that the number of electrons ne+n^ = ntot is conserved. 
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Figure 2. The LUKE momentum space is divided into two separate populations: 
the bulk electrons with momentum p < pre and the runaway electrons p > pre- The 
knock-on collisions between the populations m and m can lead to secondary runaway 
electrons. Electrons that escape the domain p < p^e by diffusion through contribute 
to the runaway electron population 


2.2. Runaway electron growth rate 


The runaway electron dynamics implemented in LUKE captures the combined effect 
of Dreicer and knock-on processes. The evolution of the runaway electron population 
under the influence of a constant electric field is calculated. Figure]^ shows the evolution 
of a runaway electron fraction with and without knock-on collisions. At first, there are 
very few runaway electrons, the knock-on contribution becomes significant only when 
a primary runaway electron population has been built up by the Dreicer effect. Then, 
an exponential growth of the runaway electron population is observed - describing the 
avalanche effect - and quickly becomes dominant over the Dreicer generation. Both the 
Dreicer and avalanche mechanisms are proportional to the bulk density = ntot — Ur, 
such that the runaway production rate can be expressed in the generic form 


dur 


Ue (7d + 7a) -t 


1 dUr 
{ntot - nr) dt 


Id + lA- 


(9) 


To quantify the avalanche growth rate, the avalanche term may be expressed as 
7 a = nr^Ai where 7 a is an avalanche multiplication factor. Thus, Eq. [^becomes: 


1 dnr 
{ntot - nr) dt 


Id + nr'^A- 


( 10 ) 


Equation 10 is an affine function of nr{t), where the constant term is the Dreicer 
growth rate and the avalanche multiplication factor is given by the slope. In Fig. 
the growth rate given by Eq. [^is illustrated for E/Ec = 40 and 60 and Tg = 0.5 keV. 
The growth rates from the LUKE calculations are evaluated numerically, the Dreicer 
as a constant value {'^d) and the avalanche multiplication factor (7 a) from the slope of 
the curve. The Dreicer growth rate calculated by LUKE agrees well with predictions 
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Figure 3. The fraction of runaway electrons {E/Ef. = 30 and Tg = 0.5 keV) as a 
function of time normalized to thermal collision time, with and without the avalanche 
effect due to the knock-on collisions. . 



Figure 4. The growth rate in constant electric field and K = 0.5 keV for E/E^ = 40 
(the curves with lower growth rate) and E/E^ = 60 (curves with higher growth rate) 
as a function of the runaway electron density, with and without the avalanche effect. 
The Dreicer contribution is in good agreement with Kulsrud’s theory m- The growth 
rates are normalized to the thermal collision frequency {vth = ^/r{vth)) 


from Kulsrud (Ref. jT 8 ]) where the Fokker-Planck equation is solved numerically. The 
avalanche multiplication factor 7 ^ characterizes the tendency of a runaway avalanche 
to develop, for a given magnetic equilibrium and parallel electric held. The actual 
runaway production due to avalanche is however time dependent since it is a product of 
the avalanche multiplication factor ja and the time dependent runaway electron density 
Ur- For example, ja can be non-zero, even though the number of runaway electrons born 
from knock-on collisions is negligible until a seed of primary electrons is established. 
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Figure 5. Avalanche multiplication factor as a function of the lower momentum cut 
off Pmin/pth for Tg = 5 keV, normalized to the avalanche factor at Pmin = Pc- 


For the sake of simplicity, the Rosenbluth model in Ref. |3] considers only secondary 
electrons born with a momentum larger than p^. However, electrons accelerated via 
a knock-on collision to an intermediate momentum pth < p < Pc could contribute 
to the runaway growth rate indirectly by populating the suprathermal region and 
thereby modifying the Dreicer flux. Numerically, three thresholds must be dehned 
when implementing the Rosenbluth model Q: the minimum and maximum values for 
the secondary electron momentum, and the minimum value Pre above which runaways 
are counted as primary electrons in the knock-on process. In order to determine these 
parameters, the lower threshold above which knock-on collisions are included is varied 
and the results are shown in Fig. [^for electric held E/Ec = 2 and E/Ec = 5 (Tg = 5 
keV). We can see that the indirect contribution of knock-on collisions to suprathermal 
energies pth < p < Pc is negligible, such that it is appropriate to set the lower threshold 
for secondary electron momentum at Pc- Energy conservation imposes that the higher 
threshold for secondary electron momentum is lower than p^e- We see that setting 
Pre = 4pc is sufficient to account for more than 80% of knock-on collisions while ensuring 
energy conservation. 

2.3. Bounce-averaged knock-on source term 

Since knock-on accelerated electrons emerge with high perpendicular momentum |3], 
it is necessary to properly account for the guiding-center dynamics in non-uniform 
magnetic held geometry and treat the full 2-D momentum electron dynamics. In a 
non-uniform magnetic held, the magnetic moment is an adiabatic invariant such that 
the guiding center parallel velocity varies along the electron trajectory. The pitch angle 
coordinate ^ in Eq. can be expressed as a function of where is the pitch 

angle measured at the poloidal position 6q of the minimum magnetic held B^igp) on a 
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magnetic flux surface. When the collisional time is longer than the bounce period [2], the 
rapid poloidal motion ensures that the electron distribution is independent 

of the poloidal angle 6. The poloidal angle dependence can thus be averaged out of the 
kinetic equation by bounce-averaging, dehned as 


Xq 


-E 
2 ^ 


d9 1 r B ^0 
27r 1-0 . f| Rp Bp ^ 


(11) 


T ^77_ 

where Rp is the major radius, and O^ax are the poloidal turning points for the 
trapped electrons. Bp is the poloidal component of the magnetic held B and the sum 
over cr applies to trapped particles (T) only. The normalized bounce time is 

1 de I r ^0 B 


with 


A(V') = 


qW = 


q{i)) 


2n 


1^ . r\Rp ^ Bp 


dO 


B 


2tt \'ijj. f\ Rp Bp 

In the code LUKE, the electron distribution is normalized to a reference density 
77,1 and the time evolution is normalized to the reference thermal collision frequency 
^coii ~ l/n(u/i), so that the resulting source term is S = S/Sl where S is from Eq. 
and Sl = Momentum is normalized to the thermal momentum p = p/pth- The 

knock-on source term is decomposed as S(p, p), 9) = — ^*{p)) where 

1 Pll - - 1 

rUeUr — 


s* = 


dvr In At *^^7(7 —1)2’ ^ ^ 

is independent of 9, so that {S} = S*{5 (.^ — ^*(p))} where f is the pitch angle cosine 
at the poloidal angle position 9 


eK„,V>.9) = 


(13) 


Here \l/('0,6') = i?('0, 6')/Ho('0) and a = sign{v\\) = sign{^o) indicates the direction of 
the electrons along the magnetic held line. Using the general relation for Dirac’s delta 
function 6(g(x)) = J2k ^{^~^k)/\g'{xk)\ where Xk are the zeros of the function g{x) and 
g\x) = dg/dx provided that g{x) is a continuously diherentiable function and g'{x) is 
non-zero: 

no-etm 


Hi-C) = Y. 




(14) 


where 91 is the poloidal angle at which the secondary electron emerges. From Eq. 13 91 
is given by 


a 




r = o, 


or 


Bg* 1 - £*2 


(1-S)(7 + 1) 


(16) 


( 16 ) 
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Using Eq. [1^ the delta function can be expressed as 

^ Bei eo 

^ \il) . 'r\ 
k 




ir 


AgTT ^ \ fi)-f\ei RpBp,e*Jei 0^)|(1 - ^o)' 

and since Be* = (1 - (1 - ^o)^o with ^e* = C 




1 rg* Bo 


1^0 


(1 _ ^* 2 ) 


A^TT ^ |-0 . f\g* Rp Bpfi 

and the normalized, bounce-averaged avalanche operator becomes 

1 1 




2 ', 2 ’ 


{S(p,' 0 ,^o)} = 


2Tr'^ In At K 


-UpUr 
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Bo 1^01 


p^'ii'y - 1) Ag (1 - ^ L|^. f I Ep |d^' 


E 


using the relation = ( 7 ^ — 1 ) = (7 — 1)(7 -|- 1 ) and p = p(3th- 


(17) 


(18) 


3. Effect of toroidicity 

The reduction of the Dreicer runaway rate away from the magnetic axis has been 
identihed in previous work jT9] including with the code LUKE, which solves the bounce- 
averaged kinetic equation [2]. At least three effects contribute to reduce the growth 
rate : the overall effect of the electric field on trapped electrons cancels out over one 
bounce period; the acceleration of passing electrons is also reduced as their pitch angle 
increases towards the high field side; the existence of a magnetic trapping cone creates 
larger pitch-angle gradients in the circulating region, thereby increasing the effect of 
pitch-angle scattering. 

As discussed in Sec. secondary electrons emerging from the knock-on collisions 
are typically highly magnetized. Since the trapping effect increases off the magnetic 
axis in a non-uniform magnetic field conhguration, the further away from the magnetic 
axis the electrons appear, the more they tend to be born trapped j3]. 

To quantify the tendency of magnetic trapping, the evolution of the runaway 
population is calculated in a scenario with circular plasma cross section and magnetic 
non-uniformity, with inverse aspect ratio ranging from e = 0toe = a/i? = l. The 
inverse aspect ratio of the Tore Supra tokamak is e ~ 0.3. The calculations in Fig. 
reveal that the runaway electron population grows significantly slower off the magnetic 
axis than in the center. 

In order to study the trapping effects on the runaway population, the Dreicer growth 
rate 7 /) and the avalanche multiplication factor 7 ^ are calculated with the bounce- 
averaged code LUKE and quantified separately. The Dreicer growth rate is found to 
be strongly affected by the non-uniformity of the magnetic field, as shown in Fig. 

A fit of the numerical results gives an analytic expression of the Dreicer growth rate 
lD/lD,cyi = 1 — 7 / 2 e/(1 -|- e). The results indicate that for e > 0.5 runaway generation 
from Dreicer acceleration vanishes. 
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Figure 6. The evolution of the runaway electron population, including the avalanche 
effect owing to knock-on collisions, depends strongly on the radial position in a non- 
uniform magnetic field configuration, where e = r/R is the inverse aspect ratio 
coordinate. E/Ec = 40, K = 0-5 keV and the time t is normalized to the thermal 
collision time Tth- 



Figure 7. Radial dependence of Dreicer growth rate, normalized to the growth rate 
for cylindrical case e = 0 and a fit yn/lD^cyi = 1 ~ l-2-\/2e/(l -|- e). 


A reduction of 7 ^ away from the magnetic axis is observed in Fig. with an 
avalanche multiplication factor that decreases with the inverse aspect ratio. In order 
to derive an analytical estimate for the avalanche growth rate including the effect of 
magnetic trapping owing to a non uniform magnetic conhguration, it is assumed that 
all electrons with momentum p > Pc will contribute to the runaway population (as in 
Ref. 0 ), except the secondary electrons that appear in the trapped momentum region 
p < Pt- The magnetic trapping criterion on the momentum pt of secondary electrons 
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Figure 8. Radial dependence of the avalanche multiplication factor from bounce- 
averaged LUKE calculations (circles), normalized to to the avalanche multiplication 
factor for the cylindrical case e = 0. The numerical integration over the knock-on source 
term in Eq. |^with the toroidal dependence in the momentum integration boundary is 
plotted in squares. The solid line shows the analytic estimate of the growth rate off 
the magnetic axis from Eq. |A.4| 


born via knock on collisions is 


m 

Brnn/r 


> 


\/l + Pr + 1 


(19) 


where Bmax/B{9) = (1 + e cos 6*)/(l — e) in a circular concentric magnetic configuration. 
Electrons will run away if their momentum exceeds both the critical momentum and the 
trapping condition in Eq. The lower integration limit Pmin for the analytical estimate 
of the avalanche growth rate is thus given by max(pc,PT)- An analytical expression for 
the inverse aspect ratio dependent avalanche growth rate is obtained by integrating the 
source term from over momentum space from Pmin to Pmax = oo , which results in a 
locally modified, inverse aspect ratio dependent avalanche growth rate 




1 1 


2 In At 
1 1 


UeUr 
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, . i^enr min —, , 

2 In At \ Ec 2e{l-\-COs9) 


( 20 ) 


The flux surface averaged growth rate is derived in Appendix Appendix A For 
eE/Ec ^ 1, 9b ^ TT , the growth rate is reduced by a factor (1 — e)^/ eE/E^j ■ 
The inverse aspect ratio dependence of the estimated avalanche growth rate obtained 


from Eq. is compared to numerical results. In addition, a numerical integration 
of the source term is also performed, with the same criterion on the lower integration 
boundary in momentum space Pmin as the analytic estimate in Eq. A.4 . The analytic 
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result is also compared to avalanche growth rate from Fokker-Planck calculations with 
the LUKE code. In that case, the trapping conditions are the same as in the analytic 
result, except for that the critical momentum is pitch angle dependent p^ = Ed (E^). 
The LUKE calculated avalanche multiplication factor and the analytical estimate show 
good agreement (Fig. [^. 

Figure|^shows the reduced growth rate for E/Ec = 5, relative to a cylindric plasma, 
equivalent to the growth rate on the magnetic axis (e = 0). Numerical integration of 


the source term shows good agreement with the analytic result (Eq. A.4). Close to 
the center, at low inverse aspect ratio, the effect of trapping is not visible, since the 
critical momentum is higher than the trapped momentum over the whole flux surface. 
This effect decreases with increasing E/Ec as the critical momentum pc decreases and 
becomes less restrictive compared to the trapping condition px-, which explains the 
flat top seen in Fig. However, for the FP calculations the magnetic trapping effect 
influences the growth rate also close to the magnetic axis. A possible explanation is 
pitch angle collisions that couple the dynamics between the trapped and the passing 
region. 

The growth rate obtained from bounce-averaged calculations suggest that the 
formation of runaway electrons is slower the further away from the magnetic axis they 
appear. In other words, the time scale of the local growth rate could be longer than 
suggested by collisional theory mm- Potential loss mechanisms, such as transport of 
fast electrons due to magnetic held perturbations jJUj could therefore act more efficiently 
on the runaway electrons formed off the magnetic axis than the ones formed on axis 
which could lead to well conhned runaway electrons at the center of the plasma. 


4. The relative importance of the avalanche effect 


The results presented in Sec. 2.2 (see Fig. have shown that the 


runaway 

electron distribution can be signihcantly modihed by including the effect of knock-on 
collisions. In order to understand the mechanisms that govern the runaway electron 
generation processes a parametric study is performed with the aim to investigate which 
runaway formation process, Dreicer or avalanche, dominates in non-disruptive tokamak 
experiments. 

The relative importance of the avalanche mechanism to the Dreicer mechanism can 
be estimated by comparing the analytic avalanche growth rate in Eq. and the Dreicer 
generation that is derived in Ref. [6]: 
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Figure 9. The fraction of runaways originating from knock-on collisions {uAlrir) 
as modelled in LUKE. The analytic estimate of when 5% (cyan line), 50% (yellow 


line) and 90% (red line) of the runaways come from avalanche is obtained from Eq. 21 


Relation to non-disruptive scenarios where runaway electrons were generated in several 
tokamaks. All data points are taken from Ref. m except for the Tore Supra (TS) 
point (discharge #40719) and COMPASS points (discharge #8555 and #8630). 


By letting a small fraction of electrons run away in LUKE, the relative importance 
of the avalanche effect as a function of plasma temperature and electric field strength 
can be evaluated numerically from the fraction of the runaway electrons that originate 
from Dreicer and knock-on collisions. In Fig. |^the fraction of runaway electrons born 
from knock-on collisions is shown, when 1% of the initial electron population has 
run away in a cylindrical magnetic configuration with constant electric field, density 
and temperature. The fraction of runaway electrons has to be small enough for the 
equilibrium parameters to remain constant. The relative importance of secondary 
runaway electrons grows at lower temperature and electric field, as the slower primary 
generation in high collisionality (low temperature) allows for runaway avalanches to take 
off. The time required to reach the runaway fraction varies strongly in the parameter 
space presented in Fig. The time required for 1% of the electrons to run away is 


illustrated for various electron temperatures (Tg = 0.05, 0.5, 2 and 5keV) in Fig. 10 


The formation of runaway electrons slows down as the collisionality increases at lower 
bulk temperature. 

The numerical results are compared with the analytical estimate from Eq. 21 with 
nr/Uf, = 0.01. The condition for the dominance of the avalanche effect '^a/id > 1 is 
plotted in Fig. along with the boundaries for which UA/rir = 5% and 90%. 

In order to relate the study to experimental tokamak scenarios, it must be 
noted that the simulations are performed for constant electric field and temperature. 
Consequently, the study is restricted to non-disruptive scenarios with well-diagnosed and 
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n=0.01n 
r tot 



E/E 


Figure 10. The time required for 1% of the Maxwellian electrons to run away, for the 
electron temperature Tg = 0.5 keV (dashed line), 2 keV (solid line and squares) and 
5keV (solid line and circles). 


quiescent conditions from several tokamaks, where runaway electrons have been observed 
in the current flattop with the relevant plasma parameters maintained essentially 
constant. Results from scenarios with reproducible measurements of electron density, 
loop voltage and plasma temperature at the runaway electron onset from DIII-D, FTU, 
TEXTOR, Alcator C-Mod and KSTAR were recently published in Ref. jH]. From 
this study the threshold electric field normalized to the critical field is found to be 
significantly higher than predicted by collisional theory where the birth of runaway 
electrons is predicted at E/Ec > 1, provided that no additional runaway electron loss 
mechanisms are present [18]. However, the condition for runaway onset in collisional 
theory does not take the time required to generate runaway electrons into account. 
Estimations from LUKE calculations in Fig. lU shows that this time scale can be 
unrealistically large as compared to the tokamak discharge duration. The time to 
generate a small fraction of runaway electrons from a Maxwellian distribution is finite 
for E/Ec> 1 but as E/Ec —)■ 1, the required time to generate runaway electrons t ^ oo. 
However, it is not the only explanation since the study in Ref. |14| found that the E/Ec 
threshold for suppression is also well above unity. 

Runaway electrons have been generated in the Tore Supra (TS) tokamak in low 
density discharges (ue < 10^® m~^). The TS pulses ^40719 and ^40721 are performed 
after a boronization and suprathermal electrons are observed in the former discharge 
but not in the latter. Both are ohmic discharges at Ip = 0.6 MA in the current 
flattop. Suprathermal electrons are observed in by the ECE edge chords at 

current ramp-up and ramp-down, when the density is low (< Uf. >= 0.4 ■ 10^® m~^), see 
Fig. |ll(a) The uniform E-field, estimated as the time derivative of the resistive flux 


21], is E» = 0.038 ±0.003 V/m and the core temperature is 3.8 keV. The determination 
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TS#40719 TS#40721 




Figure 11. Signature of suprathermal electrons on the edge ECE chord at around 
t = 14.5 s are seen in the Tore Supra discharge 47019 (left). In a following discharge, 
with higher density (right), there is no sign of suprathermal electrons. 


of the magnetic flux at the plasma boundary is described in Ref. j22]. No suprathermal 
electrons are detected by electron cyclotron emission (ECE) in the following pulse 
^40721 at a higher electron density, see Fig. ll(b)[ Similar result is found from HXR 
measurements from the vertical camera detecting emission of 20 — 200 keV (Fig. 12). A 
peak of photo-neutrons is observed at the plasma termination for the lower density shot 
{ff4Q711f) but not for the higher density shot (^40721). From the combined observations 
on ECE, HXR and photo-neutron measurements, the presence of relativistic electrons 
during the ramp-down of ^40719 is identified. During the current flattop of ^40719, 
the electron density is < n^. >= 0.64 ■ 10^® corresponding to E/Ec ~ 8, but 

there is no sign of suprathermals until E/Ec ~ 11. The suprathermal generation in 
#40719 is added to the {E/Ec^Tf) scan (see Fig. and lands in the region where 
Dreicer generation is dominant. In the higher density pulse (#40721) E/Ec ~ 4 during 
the current flattop and no suprathermal electrons are detected. These results are in 
line with those of Ref. ra where E/Ec ~ 3 — 12 is required to generate a detectable 
population of runaway electrons in the various tokamaks. 

Relating the data from the experiments in Ref. [H] and the TS discharge #40719 
to the parameter scan done in LUKE (Fig. reveals that the scenarios fall in or 
close to the region where the avalanche mechanism becomes significant for the runaway 
electron growth rate (Fig. [^. Data from two COMPASS discharges where runaways 
were observed (#8555 and #8630) fall in the region where the Dreicer effect is dominant 
|23) . Runaway electrons are commonly produced in the current ramp-up phase of the 
COMPASS tokamak, due to the relatively high E/Ec ratio (20 — 200). The circular 130 
kA discharge #8555 was part of the electron density < Ug > scan from 1 


• 4 • 10^® m 


where < Ue > for this particular shot was 2 ■ 10^^ during the flattop. The rise in 
runaway activity was observed with HXR Nal(Tl) scintillator and photoneutron detector 
as the < He > decreased from discharge to discharge, while Parail-Pogutse instability 
appeared for all discharges with < Ug > lower than in the discharge #8555. D-shaped 
160 kA discharge #8630 was done for the purpose of the sawteeth-runaway correlation 
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Figure 12. HXR data from the vertical camera (channels 1-21) in the energy range 
Ehxr = 20 — 200 keV. The HXR emission produced in the current ramp down in 
40719 is a signature of suprathermal electrons, whereas in the higher density discharge 
40721 no HXR emission is detected. 

studies with the electron density < Ue >= 9 ■ 10^® m~^. Even though the discharge had 
relatively high < >, the runaway activity correlated with the sawteeth instability 

was visible in HXR and photoneutron signals. These two COMPASS discharges ^8555 
and ^8630 are plotted on Fig. where E/Ec at the ramp-up phase were 85 and 94, 
respectively. The electron density at the time of the runaway detection is rie = 1.1 and 
0.80 m~^. In COMPASS, interferometry is used for the line averaged electron density 
< He > measurements, while Thomson scattering is used for electron temperature Tg 
and electron density Ue profile measurements. These observations suggest that knock-on 
collisions may contribute to the formation of runaway electron generation in tokamak 
plasmas, even in non-disruptive scenarios. The study predicts that avalanches can play 
an important role during current flattop. A self consistent electric field and equilibrium 
solver would be necessary to study avalanches with LUKE in disruptions, but is beyond 
the scope of the current work. 

5. Conclusion 

In this work the growth of runaway electron populations through the combined effect 
of Dreicer and knock-on collision mechanisms is studied. The Rosenbluth [3] model is 
extended to non-uniform magnetic field configurations and implemented as a bounce- 
averaged conservative source/sink term within the kinetic equation in the 3-D Fokker- 
Planck solver LUKE. Dependencies of key parameters such as electric field strength, 
electron temperature, and density are investigated. In addition, magnetic trapping 
effects are quantified in a non-uniform magnetic equilibrium, resulting in a reduced 
runaway population off the magnetic axis for both the Dreicer and the avalanche 
mechanism. 

The kinetic modelling of the formation of runaway electrons is restricted to non- 
disruptive scenarios as found in the current flattop with non-transient electric field 
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and plasma temperature. Modelling the rapidly varying temperature and electric 
field found in disruptions would require a proper description of the thermal quench 
with implemented radiative or convective loss mechanisms of the plasma energy, 
including MHD instabilities. The extension of the kinetic code LUKE to 3-D magnetic 
configuration and its coupling with a fluid code such as JOREK na would be necessary 
for such a purpose and is beyond the scope of this work. In the present paper, runaways 
are confined to the flux-surface where they are generated, such that the growth rate 
derived herein should be considered as upper estimates. 

Since knock-on accelerated electrons emerge with high perpendicular momentum, 
the full 2-D guiding-center momentum dynamics is taken into account via a bounce- 
averaged description. The effect of magnetic trapping of the electrons in a non-uniform 
magnetic field configuration, known as the magnetic mirror effect, has been investigated, 
revealing reduction of both Dreicer and avalanche mechanisms off the magnetic axis. 
An analytical expression for avalanche growth rate accounting for magnetic trapping 
is derived. It is in agreement with numerical simulations and shows that a significant 
proportion of secondary electrons are knocked into the trapped region off the magnetic 
axis. The reduction of the off axis growth rate implies that the time scale of runaway 
formation is longer at the edge than close to the center, which means that potential 
loss mechanisms such as radial electron transport could compete with the acceleration 
of runaway electrons at the edge. 

Moreover, quantifying the relative importance of avalanche generation as a function 
of plasma temperature and electric field strength, the simulations reveal that runaway 
electrons originating from knock-on collisions dominate at low temperature and electric 
field strength and likely play a significant role in runaway generation processes in several 
tokamaks with data from non-disruptive scenarios that are presented in Ref. |14( . The 
onset of runaway electrons found in these experiments is related to LUKE simulations of 
corresponding electric field and temperature in order to evaluate the importance of the 
avalanche effect, revealing that knock-on collisions may play a significant role also in non- 
disruptive scenarios. The LUKE calculations predict runaway electron generation also 
in a near critical field, in agreement with collisional theory if no other runaway electron 
loss mechanisms than collisional damping are present. However, the time to generate 
runaway electrons can be significantly large compared to the duration of the phase in 
which E/Ec > 1 in experiments. In addition, the required time for runaway electron 
formation is very sensitive to the bulk electron temperature. The lack of runaway 
electron signatures near the critical electric field could therefore be explained by the 
long time scale required for their formation. To understand this discrepancy between 
observations and theory, the existence of additional loss mechanisms that dominate 
during the current flattop must be addressed. One possible candidate is transport of 
fast electrons due to magnetic field perturbations 120]. Once such additional runaway 
electron loss mechanisms have been identified, the LUKE code may form an excellent 
test bed for quantifying these effects, which will be the objective of future work. 
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Derivation of toroidicity dependent avalanche growth rate 

As described in Sec. the avalanche growth rate is evaluated by the flux surface averaged 
knock-on source term in Eq. where the lower integration boundary is set by the 
maximum of the critical momentum and the momentum deflning the boundary of a 


passing and a trapped electron p^, given by the trapping condition in Eq. 19 For flnite 


E/Ec, the critical momentum pc> 0 . As the growth rate is averaged over the poloidal 
angle, Pmin —t Pc as the high fleld side is approached {pt —)■ 0 as 6 * —)■ tt). The growth 
rate becomes: 
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The poloidal angle Ohound where Pc = Pt constitutes the boundary between the region 
where the avalanche rate is limited either by the drag force or by the magnetic trapping 
effect. This angle is obtained from the condition Pc = Pt'- 
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If eE /{Ecil — e)^) < 1/4, pc is the lower integration limit Pmin and if eE /{Ecil — e)^) > 
1/4, Pmin = Pt{0). Averaged over the flux surface according to volumic flux surface 
average the growth rate is: 
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In the above calculation circular concentric flux surfaces are considered so that 7^1=1, 
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For eE/Ec ^ 1, (9^, ^ tt and the growth rate is reduced by a factor (1 —e)^/ (^^eE/E, 



